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We investigate the effect of a polariton bath on the photoluminescence (PL) spectrum in confined 
nanocavity-quantum dot (nC-QD) systems. We model the nC-QD system as a two-level exciton 
in strong coupling with a single photonic cavity mode interacting with its environment. The non- 
hamiltonian processes induced by the environment are taken into account via a Born-Markov master 
equation which includes: gain and loss of excitons and photons into/out of the cavity, a dephasing 
mechanism produced by phonon scattering in the semiconductor lattice and gain and loss of po- 
laritons due to the already mentioned polariton bath. In order to validate our phenomenological 
model, we fit an experimental spectrum to extract the values of all parameters appearing in the 
master equation. Our results show that polariton pumping and loss rates are comparable to the 
other parameters and therefore we have a first evidence that the polariton bath we proposed has a 
significant role in the dynamics of a nC-QD system. 

PACS numbers: 78.67.Hc, 42. 55. Ah, 42.55.Sa 



I. INTRODUCTION 

Nanocavity quantum electrodynamics (nC-QED) — the 
complex interplay between the methods of quantum op- 
tics and the specific nature of solid state systems — 
has recently emerged as one of the most active and 
promising research fields^*. A single quantum emit- 
ter, a quantum dot (QD), coupled to one confined 
light mode in a photonic crystal, micropillar, or mi- 
crodisk cavity is one of these nC-QED systems with 
more perspectives for device applications and fundamen- 
tal physics such as: single photon sources^, quantum 
information processing^, strong coupling (Sc jMH non- 
classical lighPESl, strong optical nonlinearitie d 11 * 12 ^ and 
control of cavity reflectivity^. More recently, it has been 
reported controllable all-optical switching between laser 
pulses at the single photon leve l using a QD embedded 
in a photonic crystal cavitjfl^EU. 

One of the most remarkable SC p heno mena in nC- 
QED systems is the polariton physics^^. A polariton 
quasiparticle is an entangled state between semiconduc- 
tor excitons and cavity mode photons that arise from 
the strong-coupling regime of the light-matter interac- 
tion. Despite more than twenty years of conducting the- 
oretical and experimental research, the community does 
not agree yet w ith the main physics behind polariton 
phenomen et 1 7 * 1 8 [ The large ground state occupation with 
spontaneous coherence build up can be understood in 
terms of a polariton lnsci* 1 ~* " or as a non-equilibrium 
Bose-Einstein condensat j 18 | 22 | 23 l 

Despite that SC is usually identified through their 
characteristic anticrossing in the photoluminescence 
emission spectrurrP, detailed theoretical works have 
shown that this feature is not an unique fingerprint of 
SCPl The line splitting in the spectrum depends mainly 
on the interplay of dissipative and dephasing rates and 
therefore when analyzing an experimental spectrum a 



comprehensive modeling theory is needed. A surprising 
example of agreement between theory and experiment 
was done by Laucht et al. in Ref. [2U By fitting the 
theoretical spectrum with the experimental one, they ex- 
tracted the whole set of parameters that characterize the 
interaction mechanisms between the nC-QD system and 
the environment: gain and loss of excitons and photons 
into/out of the cavity and a dephasing mechanism due 
to phonon scattering in the semiconductor lattice. 

Motivated by the polariton laser experimental 
results^, recent theoretical works have propose d 20 * 25 ! 
that some of the coherence properties of nC-QED sys- 
tems can be understood in terms of an effective pumping 
of polaritons. In Ref. [201 the authors have developed a fi- 
nite system theory of a multilevel QD which qualitatively 
reproduces the experimental polariton laser threshold. A 
Jaynes-Cummings approach was used in Ref. |55]to show 
how the incoherent polariton pumping is able to both 
sustain a large number of photons inside the cavity with 
Poisson-like statistics, and induce a separable exciton- 
photon state in the stationary limit. In addition, the 
authors also have shown that the polariton pumping is 
unable to modify the dynamical regimes of the system. 
Despite the huge literature on polaritons in nC-QED sys- 
tems, a deeper discussion of a comprehensive polariton 
bath theory is required. 

In this work, we propose a master equation with a 
polariton bath, in addi tion to the already reported pro- 
cesses in the literatur e! 20 * 24 * ^. Some of the advantages of 
this procedure are: the recognition of polariton loss as a 
relevant scattering process, the improvement in the defi- 
nition of polaritonic transition operators as discussed in 
section [II A| and the possibility to obtain a better agree- 
ment when fitting an experimental spectrum. In view of 
the high degree of complexity encountered in solid state 
systems, this could represent another step towards a bet- 
ter understanding of the nC-QD system. 
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The rest of the paper is organized as follows. In Sec. 

we describe the theoretical model and write down the 
master equation with all the processes involved. We also 
outline the Quantum Regression Theorem (QRT) equa- 
tions, which permit us to calculate the emission spectrum 
beyond the linear regime. In Sec. |III[ numerical results 
for the luminescence spectrum are shown. We investigate 
the influence of our definition of the polariton transition 
operators and discuss the different approximations to the 
first excitation manifold found in the literature. In addi- 
tion, the polariton pumping and loss rates by fitting an 
experimental spectrum are found. Concluding remarks 
are presented in Sec. [W] Finally, the appendices give 
some details about the density operator matrix elements, 
the QRT equations, and the approximations to the linear 
regime. 



represent quasi-particle states, where \Gn), \Xn — 1) are 
the "bare" states of the nC-QD system and A = lu x — ui c 
is the detuning. These quasi-particles are called polari- 
tons and their existence becomes evident by looking at 
the level anti-crossing in the emission spectr£p2l. 



A. Master equation 

In addition to loss and gain of excitons and photons, 
and a d ephas ing mechanism already considered in the 
literatur e ! 24 1 26 1 , we propose a polariton bath that accounts 
for gain and loss of polaritons into and out of the system. 
Theses processes are phenomenologically modeled by the 
system-environment interaction hamiltonian 



II. THEORETICAL MODEL 



We consider a single quantum dot in interaction with 
one photonic mode of a semiconductor nC. For the pur- 
poses of our work we treat the QD as a two-level sys- 
tem (monoexcitonic regime). This model is customary 
in the literature and captures the main physics of the 
systenPSHI] Hence, an appropriate way of describing 
the intrinsic dynamics of the nC-QD system is by means 
of the Jaynes-Cummings hamiltonian 2 ^! 



H, 



a z + hio c a) a + hg(a + a + U-a) 



(1) 



where er + , <t_ and a z are pseudo-spin operators for the 
QD with ground state \G) and excited (exciton) state 
| A") separated by the exciton energy hui x . (a) is the 
creation (annihilation) operator of the nC mode with fre- 
quency lu c , and g is the coupling strength between the 
photonic mode and the exciton. We are interested in the 
strong coupling regime. In this regime, the "dressed" 
states of the system 



com — J \Xn-l)+sm^j \Gn), 

|n, -) = - sin (^fj \Xn 1) + cos \Gn), (2) 

4> n = tan _1 (2gv/n/A), 



hS p £ = EE fi 4(i,rf+w^ )t ), (3) 



where bij n> u, b\j n R are bosonic operators for the bath 
(n) p (n)t 

fined by 



and , Pij l> * are polaritonic transition operators de- 



P. 



(n) 



,\n-l,j){n,i\ (A„->-A n _i); 



=y /a\n,i){ n -l,j\ (A n _!^A n ), 



(4) 



where i,j € {+,—}, n = 1,2, ...,oo, and A„ = 
{\Gn), \Xn — 1)} is the nth excitation manifold. Hp 
takes into account situations in which a boson from the 
reservoir is absorbed {bij n ^) causing a polaritonic tran- 
sition |n — —> \n,i) (Pij) along with the opposite 
processes. Both channels are c ontrol led by a coupling 
constant /if*„. In previous worka^ES similar transition 
operators to those on Q were defined. However, we have 
introduced the factor yjn to account for the non-linear 
behavior of polariton gain and loss rates with excitation 
number as can be expected if we look at their photonic 
analogs (see Eqs. (Al)). Additionally, the polariton 

n|n, ±), can be 
which somehow 



number operator defined by Np\n,± 



written as N P = (1/2) Pj?*P$ 

justifies our election of the -y/n factor. 

Following the Born-Markov formalis m 1 30 * 31 ! and using 
a system-environment interaction hamiltonian contain- 
ing we arrive to a master equation for the density 
operator of the nC-QD system: 
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~ = - Uh s ,p\ + y (cr z/ 9<7 2 -p) + -y ( 2at P a - aa) p- paa)) + y (2apa) - a) ap - pa) a) 

Pa 7a 
+ -^{2(T+PV- - (T-V+P - PP-<T+) + y(20-_p(T + - <T + (J-p - p(T + <J-) (5) 



p„ 

ij n ijn 



+ y E - ^? )f p - pp^p^) + 1 E (api B) /4 B)t - it )f ^V - 



In this master equation, 7^ is a pure dephasing rate in- 
duced by constant scattering processes between the QD 
and the semiconductor lattice phonons, P a represents the 
photonic pumping rate caused by the emission of weakly 
coupled QD's present in the nC, 7 a is the cavity loss rate 
(coherent emission), P a is the rate at which excitons are 
incoherently pumped by optical or electrical excitation, 
7^ is the exciton decay rate (spontaneous emission), and 
finally P p and j p are respectively, the rates at which po- 
laritons are created and lost (see Fig. [IJ. 

Eq. ([5]) can be transformed in an infinite set of lin- 
ear differential equations for the density operator ma- 
trix elements on the direct product basis {\Gn}, \Xn}}. 
When the system is initially in its ground state, there 
is an independent set of equations for the nC-QD pop- 
ulations {pGn,Gn, PXn,Xn) and coherences between lev- 
els in the same excitation manifold A n+ i (pon+i.Xn)- 
This set of linear equations is found in Appendix [A| (Eqs. 
Al ). It is useful for practical computational purposes and 
it provides a better understanding of the effect caused 
by system-environment interaction processes. From Eqs. 



( Al ) it becomes clear that the overall effect of these pro- 



cesses is to cause transitions between excitation mani- 
folds and to reduce coherence in the system. 

B. Emission spectrum 

When dealing with experimental samples of the nC- 
QD system, the emission spectrum is one of the most 
accessible observables. Assuming that the principal cause 
of light emission of the system is due to cavity leakage, 
the spectral function in the stationary state is given by 



/>oo 

S(uj)(x lim Re / dre'^^ (t + tW*)), 
Jo 



(6) 



which is the Fourier transform of the cavity field cor- 
relation function G {1 \t, T ) = (a^(t + r)a{t)) weighted 
by the factor e~ 1T to account for additional broadening 
of the spectral lines due to the finite resolution of the 
detectors^. 

Our calculation of the correlation function closely re- 
sembles the one by Tejedor and co-workers^. First, we 
write operator a) as 
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FIG. 1. (Color online) Transitions between "bare" energy 
levels induced by the hamiltonian and non-hamiltonian pro- 
cesses. Continuous green arrows represent transitions due to 
the nC-QD interaction. Dot-dashed red arrows show photon 
gain and loss transitions. Dashed orange arrows depict tran- 
sitions caused by the exciton pumping and loss processes. Fi- 
nally, the dotted blue arrows correspond to the polariton bath 
induced transitions (these transitions can be better visualized 
using the "dressed" states basis). The dephasing mechanism 
is not displayed since it does not cause transitions between en- 
ergy levels, it only modifies the non-diagonal density matrix 
elements. 



i(4 



where a Gn = \Gn + l)(Gn\ and_a^„ = \Xn + l)(Xn\. Us- 



Gn 



,t 



(7) 



ing the master equation (Eq. ([5|) it is possible to find an 
independent set of linear differential equations including 



the expec tati on values of a' Gn and a Xn (see Appendix 
|B} Eqs. ( |B1[ )). According to the Quantum Regression 
Theorem (QRT), time evolution for expectation values 



(a Gn (t + T)a(t)}, (a x Jt + r)a(t)), (4 (t + r)a(t)), and 



(<T n (£ + T ) a (t)), is also governed by Eqs. (Bl ) with initial 
conditions 



( a Gn( i ) a ( < )) = Vn+lpGn+l,Gn+l(t), 
( a Xn(*) a (*)) = Vn + ip Xn +l,Xn.+l(t), 

(a-t(t)a(f)) = Vn + lpGn+i,x n {t), 

(? n (t)a(t)) = \/npXn,Gn+l(t), 



(8) 



where t stands for the time in which the system 
has reached the stationary state, and operators = 
\Xn) (Gn\ and <, n = \Gn + 1) (Xn — 1| were introduced in 
order to obtain a closed set of operators. After solving 
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can compute the first-order correlation function as 
G^(t,r) = J2V^Tl ((4(t + T)a(t)) 

n 

+ <4n(* + 1>(*)} > 



Eqs. (Bl) with initial conditions given in Eqs. we inition. For low polariton rates (P p = meV/fi), 

the difference between both spectral lines is almost neg- 
ligible. This is expected since for low pumping rates, 
the population for excitation manifolds higher than Ai 
are much less than one, making the dependence with n 
almost irrelevant. However, for higher pumping rates 
(P p = 10~ 2 meV/h), both spectral functions appreciably 
differ from each other. In fact, the interior peaks have a 
greater intensity if the spectrum is calculated with lad- 
der operators given by Eqs. Q. A gain, this is expected 
since for n-dependent polariton pumping, the second ex- 
citation manifold populations are higher and the interior 
peaks correspond to transitions |2,±) — > |1,±). Addi- 
tionally, there is a minor increment in the line widths, 
since decoherence effects are larger when the factor \fn 
is introduced (see Eqs. (All). 



(9) 



Now, with the knowledge of G^(t, r), the emission spec- 
trum is obtained via Eq. {6}. The main advantage of the 
previous approach is that we can calculate the emission 
spectrum beyond the first excitation manifold without 
doing any further approximations. 



III. NUMERICAL RESULTS 



A. Polariton pumping and emission spectrum 



In this section, we show the effect in the emission spec- 
trum caused by the excitation manifold dependent factor 
in the polariton ladder operators definition (Eqs. Q). 
In previous workgUM25J this factor was not introduced 
which then leaves the effective polariton pumping and 
loss rates independent from the excitation number. From 
Eqs. (All, it can be seen that effective photonic rates 



grow with the photon number n, a behavior that we re- 
produce for polariton rates by introducing the factor ^fn 
in Eqs. Q. 

In Fig. [2j we show emission spectra calculated with 
and without the factor y/n in the ladder operators def- 



B. Approximations to the emission spectrum 

The importance of computing the emission spectrum 
using the exact dynamics of the first-order correlation 
function governed by Eqs. (Bl ) is developed in this sec- 
tion. For this purpose, we fix the polariton pumping and 
loss rates to zero. 

There are at least two ways of finding an approxima- 
tion to the emission spectrum. The first one is to approx- 
imate the first-order correlation function up to the first 
excitation manifold^. This leaves us with the following 
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Eqs. (4) 
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FIG. 2. Emission spectra using two definitions for the 
polariton ladder operators for P p — 10~ 3 meV/h (top) and 
P p = 10 _2 meV/ft (bottom). Continuous shadowed lines rep- 
resent spectral functions calculated with Eqs. Q, whereas 
dashed lines correspond to emission spectra calculated with- 
out the factor yfn as is used in Ref. 1251 The other parame- 
ters in the master equation were fixed to u) c = 1000meV/7i, 
g — lmeV/h, 7 a = O.lmeV/ft, and the other ones to zero 
A = P a = P, = 7ct = 1<t> = 0. 
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FIG. 3. Emission spectra calculated with three differ- 
ent methods for the parameters reported by Finley and co- 
workers^ 1 ( ff = 59 fieV/h, P a = 4.5 fieV/h, P a = 0.5 fieV/h, 
7„ = 68.0 fieV/h, 7 CT = 0.2 jteV/fi, and ^ = 19.9 jteV/ft). 
The continuous line is t he P L spectrum computed with the 
exact expansion (Eqs. (Bl|) up to A25. Dashed and doted 

1 101 



lines correspond to the spectrum calculated with Eqs. 
and Eqs. ( |TTj ), respectively. 
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set of differential equations 

^ r (ot(r)) = (i Wc -r 1 )( a t(r))+v< - + (r)) 
f (10) 
-(a + (r)) = t S ( a t( T )) + (*w a - r'O (a + (r)), 

where Ti = 3P a /2 + 7a /2 + P CT and r' x = P a + P CT /2 + 
7a-/2 + 70. Then, the first-order correlation function 
G^'(t, t) can be found by using the Quantum Regression 
Theorem with the appropriate initial con dition s. The 
second way is due to Finley and co-worker o I I , and it 
combines the first approximation with the exact solution 
as we show in Appendix [Cj Again, the problem becomes 
two-dimensional and the spectrum is computed via 

jVM) = (iw c - T 2 ) ( a t(r)) + i 9 (a + (r)) 
f (11) 
-(a + (r)) = ig{a\r)) + ^ x - T' 2 ) (a + (r)), 

where T 2 = j a /2-P a /2 and T' 2 = P a /2+ la /2+ 70 . Both 
approximations are good for low pumping rates. How- 
ever, as displayed in Fig. [3j they are unable to reproduce 
the exact spectrum for typical parameters controlling the 
dynamics of the nC-QD system in real experimental se- 
tups. 

C. Theoretically fitting an experimental emission 
spectrum 

To validate our model of the polariton bath, we fol- 
lowed a similar proc edure to the one developed by Finley 
and co-worker d l 33 l By fitting one of their experimen- 
tal spectra^ we were able to find values for the coupling 
constant g and the gain and loss rates involved in our 
master equation (Eq. Both the detuning and the 

cavity mode were experimentally determined by Finley's 
Group so they are not included in the fitting process and 
remain fixed at A = —50 fieV/h and w c — 1207.1 meV/h. 
We also use their estimated value of the broadening pa- 
rameter appearing in Eq. (rob, 7 = 30^eV/?P^. 

Based on the results of the previous section, when fit- 
ting an experimental emission spectrum it is of great im- 
portance to compute G t ^ 1 \t 1 r) using the exact expansion 
(Eq. ([9])) since for some range of parameters the two ap- 
proximations mentioned above are not able to reproduce 
the actual spectral lines. As a consequence, we do not 
have an analytical expression for S(ui) so we cannot use 
a regular least-squares algorithm to fit the spectrum. To 
overcome this obstacle we used a Simulated Annealing 
(SA) algorithm 34 defining the cost function as the sum 
of square differences between the experimental and the- 
oretical spectral lines 

S 2 = 5^(5(wi) - 5 eap (cJi)) 2 , (12) 

i 

where S exp (uji) is the experimental PL intensity at fre- 
quency Wi, and the index i runs over the available data. 
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FIG. 4. (color online) Theoretical fit and control spectra. The 
continuous red line shows the theoretical PL spectrum fitting 
the experimental data (black dots) at A = — 50/ieV/ft. The 
dashed yellow line is the control spectrum and the gray tri- 
angles are the experimental intensities at A = +50 fieV /h. 
The cost function is 5 = 12.6cts/s in the first case, and 
5 = 29.8cts/s in the second one. The fundamental mode 
frequency is fixed at uj c — 1207.1 meV/ft. 

SA has proven to be a fairly good tool to optimize func- 
tions in high parametric spaces making it the best choice 
for our purposes. 

The best fit (displayed in Fig. [4]) was obtained for g = 
67.93 (j,eV/h, P p = 0.02 fieV/h, P a = 2.27 fieV/h, P a = 
2.32/xeV/ft, 7 P = 8.96 fxeV/H, j a = 31.82 neV/h, j a = 
10.12 fj,eV/H, and 7^ = 7 AO /j,eV / h. These values show 
that polariton pumping and loss rates are comparable to 
the other parameters, indicating that the polariton bath 
we proposed has a significant role in the dynamics of the 
nC-QD system. A conclusion that can be drawn given 
the quality of our fit. 

Additionally, to strengthen our result we use a con- 
trol experimental spectrum from Ref. [24] with A = 
+50 fieV/h. The experimental data along the theoreti- 
cal curve computed with the same fitting parameters, are 
also illustrated in Fig. |4j This control spectrum shows 
that the parameters obtained are also good to describe 
the characteristic asymmetric double peak features near 
resonance. Note that the small blue shift in the spectral 
lines is consistent with the eigenenergies of the "dressed" 
states. 



IV. CONCLUSIONS 

To summarize, we have presented a phenomcnologi- 
cal model for a polariton bath interacting with a nC- 
QD system. By theoretically fitting a photoluminescence 
spectrum, we were able to show that our model agrees 
with the experiment for non-zero values of the polari- 
ton pumping and loss rates. These values have the same 
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order-of-magnitude of the rest master equation parame- 
ters, showing that at least from an experimental point of 
view the polariton gain and loss processes are feasible in 
real nC-QED systems. Additionally, we have shown the 
importance of computing the PL spectrum with the ex- 
act dynamics of the first-order correlation function. The 
two approximations we have discussed, fail to give an 
accurate description of the nC-QD system in real exper- 
imental conditions. 
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Appendix A: Density operator matrix elements 
equations 

Denoting matrix elements by pi m ,Jn — {I m \p\J n ) i 
with I,J £ {G, X}, and using Eq. (5), is straightfor- 
ward to show that 



-JlPGn^Gn — n (P a + Pp) PGn-\,Gn-l + n Pp PXn-2,Xn-2 



d 

-{(n + l)P a + nj a + P a + 2(n + l)P p + (2 - S n i)nj p }pG n ,Gn + igV™ PGn,x n -i 

-ig\fn pXn-l,Gn + (n + l)(7a + 7p) PGn+l,Gn+l + {icy + {n + l) 7p } PXn.Xn , (Ala) 



^PXn,Xn = {Pa + (jl + l)Pp] PGn,Gn + {nP a + (u + l)Pp} pXn-\,Xn-\ 



-{(n + l)P a + H'Ja +7 CT + 2(n + 2)P p + (2 - 5 n0 )(n + l)l p }pXn,Xn ~ igVn+ 1 pGn+l,Xn 
+igVn + lpxn,Gn+l + (n + 2)j p pGn+2,Gn+2 + {(« + l)ja + ( n + tylp} PXn+l,Xn+l , f^lp) 



j, PGn-\-l,Xn 

a/ n(n + 1) P a PGn,Xn—l 

+ y/{n + l)(n + 2) 7 a pGn+2,Xn+l 

+ {*A - (2n + 3)^ - (2n + 1)^ - y - ^ - l4> - 2(n + 2)P P - (2 - S n0 )(n + l) 7p | p Gn +i,Xn 

+igV n+1 pGn+l,Gn+l - igVn + 1 px n ,Xn , ((AlJ;) 

Since pxn,Gn+i — Pcn+i Xm the equation governing time evolution for pXn,Gn+i is found by taking the complex 
conjugate of Eq. (|Alp). As mentioned above, this linear set of differential equations is quite useful to find the time 
evolution of the density operator p. 

Appendix B: QRT equations 

Expectation values for a,Q n , a} Xn , a\, q n can be written in terms of density operator matrix elements using the 
master equation (Eq. [5]). Their time evolution can be easily found to be 



d_ 



(4nW> ={^c - ^(2n + 3) - ^-(2n + 1)-P a - P v {2n + 3) - ^[(2 - 8 nl )n + (2 - S n0 )(n + 1)] }<4 n (r)> (Bla) 
+P a ^n(n + 1) (a^-iW) - Wn<fe(r)> +y a \/(n + l)(n + 2) (4„+i(t)> +7<r(akW) + *<?V™ + 1 <<4(t)>; 



^( a kW) = [i"c - ~ (2n + 3) - ^(2n + 1) - 7< r - P^n + 5) - ^ [(2 - 5 n0 )(n + 1) + 2(n + 2)] <|B1 
+P a Vn(n + 1) (ak-iW) + wVn + 2 (c„+i(r)} + 7a V(n + l)(n + 2) <<4„ +1 (r)> + P ff <a^(r)> - ipVn + 1 <<t£ +1 (t)> 
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d_ 



<4W) = - Pa(n + 1) - 7a n - ^ - ^ - 70 - P p (2n + 3) - ^[(2 - 5 n0 )(n + 1) + (2 - S nl )n] }(<r£(r)) 



mb) 



j-(?n(T)> = |i(2w c - - P a (n + 1) - j a n ~ -y ~ y 



- 7^, - P p (2n + 3) - ^[(2 - <5 n0 )(n + 1) + (2 - <5 nl )n] | (? n (r)) 



P a \Jn 2 - 1 (? n -i(r)) +7a\/w(n + 2) (?„+i(r)) + igVn + 1 (4„-i(t)) - igVn (aL( T ))' 



IBli) 



From Eqs. (Bl), it becomes clear that the position of 
spectral peaks solely depends on j, uj x , uj c . This means 

that dissipation and pumping rates only modify their d , , _ f _ _ ja_ _ 1 V^/^f / r \\ 

intensity and width. This is expected since population dr ° + T 2 2 ^ J ^ °" 

dynamics and decoherence properties are ruled by these ^ (C2) 

rateS - -ts^V^^W) - (at G »», 

ra=0 



Appendix C: Approximations to the emission 
spectrum: Derivations 

Both approximations to the emission spectrum dis- 
cussed in section HI B| can be derived from Eqs. (Bl I, 
which describe the exact dynamics of the first-order cor- 
relation function <?W(t, t). The first approximation is 



straight forward. Up to Ai, we have 



'■GO 



) and 



Bl) 



(cr_|_) ~ (ctq) so that other expectation values in Eqs. ( 
become zero, leaving the two-dimensional system in . 
( 10 1. To derive the second approximation we need to find 
the exact time derivatives of (at) and (cr+). From Eqs. 



bqs. 



(Bl I, summing over all values of n and simplifying, we 



find 



It is worth noting that this pair of equations are exact and 
cannot be found by doing the first excitation manifold 
approximation. From Eq. ^ and since a + = Yl n a h> 
we have that Eq. (CI) is precisely the first line in Eqs. 



(11). However, to obtain the second line, it is necessary 



to approximate the second sum in Eq. (C2) to A 1; i.e., 



£ Vn+l({J Xn (T)) (al n (r))) « -<4 ), 



(C3) 



dr 



(at(r)) 



lLO c 



E 

n=0 



2 



To 
2 



V^+T((4W) + (»kW» (ci) 



This may be objectionable, since to find Eq. (C2) we 



had to sum over all excitation manifolds and now we are 
approximating the last sum to Ai. 
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